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Abstract 

We study the pion structure function in nuclear medium using the Nambu and Jona- 
Lasinio model, and its implication for the nuclear pion enhancement of the sea quark 
distribution in nuclei. By using the operator product expansion, medium effect of 
the nuclear matter is incorporated in calculations of the twist-2 operators. We find 
density dependence of the pion structure function is rather weak around the nuclear 
matter density. We also discuss how the medium modification of the pion structure 
affects the sea quark enhancement in the nucleus. 
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1 Introduction 



Modification of tlie liadron structure in nuclear medium is an interesting topic of 
the nuclear physics, and has attracted considerable interests. In the deep inelastic 
scattering, such modifications of the nucleon structure in nuclei are observed experi- 
mentally, that is, EMC effect It is very difficult to calculate the hadron structure 
function itself by using QCD directly, though the measured momentum dependence 
of the structure function is consistent with predictions of the perturbative QCD. Re- 
cently, several attempts were done to calculate the structure function in terms of low 
energy effective theories of QCD, which reproduce hadron properties such as masses 
and magnetic moments^. In those studies, it is assumed that calculation of the 
twist-2 matrix elements within the effective theory gives structure function at certain 
low energy scale fi ~ IGeV, where the effective theories make sense. The resulting 
structure function is evoluted to the experimental high momentum scale with the 
help of the perturbative QCD, and comparison with experiments can be made[0, ||. 

In particular, the pion structure function is of special interest so as to study 
the nuclear pion enhancement of the sea quark distribution in the deep inelastic 
scattering[p, ^ and the Drell-Yan (DY) process0. The nuclear pion enhancement due 
to the attractive nucleon-nucleon interaction in the nucleus was extensively studied in 
the last decade^]. The nucleon-hole and delta-hole excitations shown in Fig.l modify 
the pion-like mode in the nucleus, and lead an excess of the pion number in nuclei. 
Llewellyn Smith 0, and Ericson and Thomas]^ discussed that such an enhancement 
of the pion number changes the sea quark distribution in the nucleus. Assuming that 
a part of the nucleon sea quark distribution is originated from the virtual pion cloud, 
the sea quark distribution of the nucleus is written by the convolution formula^; 

Q^nn{x,Q^)= dy/yf^^^{y)q^{x/y,Q^) , (1) 

where q^j^^i^x^Q"^) is the sea quark distribution of the nucleon in the nucleus orig- 
inated from the pion cloud, and q-nix-niQ"^) the antiquark distribution in the pion. 
Following the work of Ericson and Thomas, the longitudinal pion momentum distri- 
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bution in the nucleus f^^Mhj) is given by 



fl'^ poo rq—mj\[y 



where Qttnn is the pion nucleon couphng constant, F^^^{t) the pion-nucleon form- 
factor, and m^v, the masses of the nucleon and pion, respectively. Here, R{uj,q) 
is the response function of the nuclear matter calculated by the random phase ap- 
proximation, and carries the information of the collective excitation illustrated in 
Fig.ip. In the case of free nucleon, we basically replace R{uj, q) with 1, and recover 
the Sullivan's convolution formula for the single nucleon. 

By using substantial enhancement of the sea quark distribution in nuclei was 
found by several authors]^, 10], and such an enhancement should be observed in the 
DY experiment 0]. However, recent experiment of the DY process indicates that there 



is no enhancement of the sea quark in nuclei|]TT]]. Also, the sea quark enhancement 
due to the pion excess is not observed in the deep inelastic scattering[Q. Recently, 
Bertsch et a/.[jl3l and Brown et al-^^ studied such a puzzle 'why the nuclear pion 
enhancement is not observed', by taking into account the modification of the gluon 
properties in the nucleus, or the universal scaling law due to the chiral symmetry 
restoration in the medium ||15||. 

In the above discussion, the medium effect leads an excess of the number of pion 
in nuclei, but the pion structure function itself is assumed to be unchanged in the 
nuclear matter. In this letter, we concentrate on the medium effect on the quark 
distribution in the pion, and its implication for the nuclear pion enhancement. If the 
medium effects change the pion distribution in the nucleus, it is also important to 
examine whether the nuclear medium modifies the quark distribution in the pion or 
not. Following this scenario, we study the density dependence of the pion structure 
function. 

The antiquark distribution of the pion appeared in (|l|) consists of the valence 

antiquark and sea quarks. For x > 0.l[|, most contribution to the convolution formula 

^ We do not consider the region x < 0.1, since the nuclear shadowing effect, which we do not 
take into account in this letter, is important in that region. 
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comes from the valence antiquark part. Hence, we consider only the valence quark 
distribution in the pion and its medium modifications in this letter. 

In the previous studies of the present author |T^, we have evaluated the pion 
structure function using the Nambu and Jona-Lasinio model (NJL) and ob- 
tained reasonable agreements with experimental data. Here, we adopt the similar 
procedure in the nuclear medium. The NJL model, where the chiral invariance is 
main ingredient, describes the dynamical mass generation of the constituent quark, 
and reproduces SU (3) meson properties successfully as a low energy effective theory 
of QCD||19|. It is believed that restoration of the chiral symmetry takes place at high 
density. Even at the nuclear matter density p = Po = 0.17fm~^, the chiral symmetry 
may be partially restored. We will show how the pion structure function changes in 
the nuclear medium in terms of the NJL model. 



2 Nambu and Jona-Lasinio model in medium 

Before we start calculating the pion structure function in medium, we briefly give 
a basic procedure of the NJL model at flnite density. The NJL model demonstrates 
the spontaneous chiral symmetry breaking of the QCD vacuum, in which the gluon 
degrees of freedom are assumed to be frozen into a chiral invariant 4-quark interaction. 
The NJL lagrangian is written as, 

Cnjl = ^Pi^r^^-m)^lJ + Gsm,^|Jy + {^p^l5t^^y] , (3) 

where -0 denotes the quark fleld with current mass m, tj are SU{3) flavor matrices 
with the normalization tr[tjtj] = 6ij/2, and Gs the coupling constant. In this model, 
the constituent quark mass M and the quark condensate (ipip) have non-zero values 
due to the dynamical chiral symmetry breaking, if the coupling constant Gs is larger 
than a certain critical value. As a result, the zero mass Goldstone pion appears in 
the chiral limit m = 0|jl8[. Using the realistic current quark masses, m„ ~ 5MeV and 
rris ~ ISOMeV, SU (3) meson properties are well reproduced in the generalized NJL 
model with the U{1)a anomaly |I9[ . 



For the application at finite density, we adopt the quark propagator; 



S{p) = {i) + M) 



1 IIT 



p2-M2 + ^£ E 



(4) 



where E = y/jP^+M^. The quark Fermi momentum kp is related with the baryon 
number density. Restricting the 2-fiavor, the baryon number density p is given by, 



r 

p = 4 

Jo 



(27r)3 • 

The quark self-energy is modified by the presence of other quarks in medium. 
From (I), the quark condensate becomes, 

(^^) = -itr J ^^S{p) 

Then, the Gap equation depends on the Fermi momentum. 

M = m-Gs{ipij) (6) 

Apparently, the constituent quark mass M vanishes in the limit kp ^ (x. In this let- 
ter, we use the 3-dimensional sharp cutoff scheme. Choice of the 3-dim. sharp cutoff 
makes it possible to maintain the gauge invariance in the following calculations [p!9[] . 

The pion properties are also obtained by solving the Bethe-Salpeter (BS) equation. 
The quark- ant iquark scattering matrix in the pion channel with the total momentum 
p is evaluated as, 

^(^^ = 1 - 2GsU{f) 

where 

d^k 



n(/) = ^tT j j^,[il,S{k)i^,S{k + p)] 



i2ny 

The pion mass is determined by the pole of denominator at p"^ = m\ 
1 - 2GsIi{p' = ml)=Q 
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Straightforward calculation yields, 
where 



Form the scattering matrix (|^), the pion-quark-quark coupling constant is ex- 
pressed as, 

^'^■^ ^ ^5Pw)] ■ 

We use the parameter set, m = 5.8MeV, A = 615MeV, G^A^ = 9.08, which are 
chosen to reproduce the pion mass 140MeV and the decay constant f^^ = 93.3MeV at 
zero density. Results for the constituent quark mass and the quark condensate are 
M = 358MeV and {tjj-ip)^^^ = — 245MeV, which are consistent with phenomenological 
values. All the results are basically same as those obtained in previous studies [1^ . 



In the nuclear medium, the quark condensate decreases as the baryon number 
density increases, and hence the constituent mass and the pion decay constant also 
decrease. At the normal nuclear matter density, po = 0.17fm~'^, values of these 
quantities become about 70% of their zero density values. The pion mass is almost 
independent of the density before the chiral phase transition. After the chiral sym- 
metry is restored, the pion mass grows rapidly. Such a behavior is general tendency 
of the NJL model calculations [pO|] . 

3 Calculation of the moments of the structure func- 
tion 

We shall consider the medium modification of the pion structure function using the 
NJL model. In order to calculate the pion structure function, we adopt a formalism 
based on the operator product expansion (OPE) to separate the hard and soft parts 
of the matrix elements explicitly ||21||. By doing so, it is easy to deal with the effects 
of finite density on the matrix elements. 



As a result of the OPE, the forward scattering amphtude is written by the hard 
part (Wilson coefficients) and the matrix elements of local operators as illustrated in 
Fig. 2. The twist-2 operators of the un-polarized structure function are defined as, 



Omi/.2-/^„ = ^V'7mi-^/^2 • • • -^MnV' + permutations (11) 

where is the covariant derivative of QCD. Here, we write the quark part only. 
Also define the matrix element A^, 

(p|CMiM2-Mn \P) = ^nPMiPM2 • • -P^n + t^acc terms (12) 

where \p) is the pion state with the momentum p. Wc also write the momentum of 
the virtual photon q. u = p ■ q, and the Bjorken-x; x = —q^/2mT^v. This matrix 
element is related with the moment of the structure function; 

dx x--'F,{x, q') = AMq']) (13) 

Jo 

where the structure function F2{x) = Y^elxqi{x) with the qi{x) being the quark 
longitudinal momentum distribution function in the pion. 

Here, we choose the axial gauge q- A^ — Hence, multiplication of q^^ - ■ ■ q^^ for 
both sides yields, 

Ar,u^ = {p\,p /l{q-dr-',p\p) (14) 

Then, we can get expressions for the n-th moment of the structure function. 

We shall evaluate the matrix element of local operators in terms of the NJL model. 
The diagram Fig. 2 indicates only the valence quark contribution up to the leading 
order. Hereafter, we stress on the n-th moment of the valence quark A^"-\ which is 
related with the valence quark distribution as, 

AZ""^ = dx x^'-^qix). (15) 

^0 

Throughout the following calculations, we omit the overall charge (isospin) factor for 
simplicity. The valence antiquark also has a similar contribution. 
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A\ is written by, 



1 /'m y.gc, y (27r)4 (A;2_M2 + Z£)[(A;+p)2-M2+^£]2 



After some algebra, we find 

= -pA,-^[p'Hp')] 

= P,- (17) 

where we have used eq. (0) for g^^qq in the last equality. Thus, we get A\°-^ = 1, which 
manifests the Adler sum rule for the valence quark distribution function. Antiquark 
contribution also gives the same result. 

We next consider the second moment which gives the momentum fraction carried 
by quarks. The 2-nd moment is defined by, 

AT^p^iPu = ^[{p\i'lfi{k + p)ui'\p) + permutation] , (18) 

Similar manipulation with the case of the 1-st moment gives, 
d'^k k^ku + k^py P,i{k + p)u 



+ 



(27r)4 ' [{k + pf - M2 + ief {{k + pf - + ie^ 
d p^ku 
^^^df (fc2 - M2 + ie){{k + pf - M2 + is) 
d p^ 
^^"^^ df (fc2 - M2 + ie){{k + pY - M2 + le)^ 

The first and second term give the trace terms, which are proportional to g^j_y, and 
are irreverent in the Bjorken limit [^]. Sum of the third and forth terms becomes, 

AT^Pi^Pv = -p^^Puglgq^p^h 



c/p2^^ J (27r)4(fc2-M2 + ie)[(A; + p)2-M2 + ie] 



-Pt.Puglqq^[p\h - \h)] 
1 



Pt^Pv ■ (19) 
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Again, the calculation for the antiquark part gives 1/2. This is just the momentum 
sum rule. 

r-l 



/ dxx[q{x) + q{x)] = 1 
Jo 



Above result shows that the momentum sum rule is strictly satisfied in the NJL 
model, i. e. all the momentum in the pion is carried by the valence quarks at the low 
energy model scale. Of course, if we consider sea quark degrees of freedom, which are 
identified with higher order loop corrections in the NJL model[0, the sea quarks 
carry the certain amount of the momentum fraction even at the model scale. 

In the previous studies |I6|, p!7| , we had adopted the covariant parton formalism and 
artificial cutoff procedure, which breaks the gauge invariance. The correct treatment 
described here maintains the momentum sum rule. 

The results for the 1-st and 2-nd moments are quite reasonable. Characteristic 
feature of the non-perturbative model may be seen in higher moments. It is worth 
mentioning that the 1-st and 2-nd moments are independent of the density. 

Calculations for the higher moments require more tedious manipulations. The 
3-rd moment is defined by. 



Pt^Pi^Pj = ^[(pI^7m(^ + P)i^if^ + P)ii'\p) + permutation] . 
We show only the result. 



+ 9l,^Mp\M) 



where 



Stt 



4(M2+A2) 



(20) 



(21) 



h{p\M) 



IK, 



4M2 



4(M2+fc2 ) 



X 



8p' 



log 



Vp^ + 



1 1 



4^2 



(22) 



Note that the 3-rd moment explicitly depends on the NJL model parameters. When- 
ever the density increases, the 1-st and 2-nd moments are unchanged, i.e. they are just 
numerical numbers. However, the 3-rd moment shows density dependence through 
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the constituent quark mass, coupling constant, and pion mass. The restoration of 
the chiral symmetry drives the reduction of the 3-rd moment, which will be shown 
later. 

Similarly, the 4-th moment is calculated as. 



< = i + 9;«S«P^A') ■ (23) 



Higher moments are calculable in the same manner. 



It is interesting to note that numerical numbers appeared in eqs. (^) and (P^), 
1/4 and 1/8, are the same as results of the delta function quark distribution q{x) = 
6{x — 1/2). Such expressions are desirable. Remaining model dependent parts are 
corrections due to the non-perturbative bound state structure. In the limit of high 
density p ^ oo, the correction terms go to zero, and the structure function approaches 
the non-relativistic (non-interacting) limit q{x) = 6{x — 1/2). 

4 Numerical results and discussions 

Using the moments obtained above and performing the inverse Mellin transfor- 
mation (or trial and error), we arrive at an expression of the quark momentum distri- 
bution function q{x). Practically, it needs calculations up to the 4 ~ 6-th moments 
to determine a shape of the structure function by using eq. (|T^). The calculated 
distribution function is defined at the low energy scale, where the effective theory 
is supposed to work. To compare with experiments, we carry out the evolution 
of the n-th moments by using the perturbative QCDpT[. The result of the effective 



theory plays a boundary condition for the QCD evolution. Here, we choose a low 
energy scale n = 0.35GeV, which is a value of the constituent quark mass, and the 
QCD scale parameter Aqcd = 0.25GeV. It is difficult to estimate ambiguities of the 
QCD evolution from such a low momentum scale. It needs more detailed studies to 
clarify use of the perturbative QCD below 1 GeV. 

We show in Fig. 3 the valence quark distribution function at 8GeV^ The solid 
■I- Of course, the valence antiquark distribution has the same shape. 
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curve denotes the model calculation at zero density; p = 0, which is to be com- 



pared with the experimental fit indicated by the dashed curve P3[|. The result is in a 
reasonable agreement with the experiment. 

We next show the results at the nuclear matter density p = pq in Fig.4 and Fig. 5. 
To deal with the medium effect, we use two different methods. 

SETl : Simple application of the NJL model at finite density. 
SET2 : Brown-Rho scaling relation 

To calculate results of SETl, we simply use the NJL model at finite density, 
described previously. Density dependence of the various properties is calculated by 
solving eqs. (H) and (|^). As the matter density increases, the quark condensate, 
constituent quark mass, and pion decay constant decrease about 30%, whereas the 
pion mass and the pion-quark-quark coupling constant are almost unchanged. Such 
a behavior was already studied by several groups, and details are found in ref. pil 



However, the confinement of quarks are absent in the NJL model, and so this 
model at finite density implicitly assumes the 'deconfined' quark matter in the nu- 
cleus. Hence, the result obtained by the NJL model at finite density is not so reli- 
able. Here, we adopt another possible method to incorporate the density dependence. 



namely, universal scaling law proposed by Brown and Rho|T^. Thus, we assume the 
following scaling relation; 



M* Q* f* 



(24) 



M g^gg U 

where the asterisk denotes the values of quantities in the medium. In this case, the 
pion mass is assumed to be constant. We use M* /M = 0.8 at the normal nuclear 
matter density. 

Results for the quark distribution function at the normal nuclear matter density 
are shown in Fig.4. Ratio of the quark distribution in the medium to one in the 
vacuum is also shown in Fig. 5. As for SETl indicated by the solid curve, density 
dependence of the quark distribution is weak. The calculated distribution function 
is almost same as one of the free pion shown in Fig. 3. Only around x ~ 1, the 
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quark distribution function decreases about 10%. Indeed, values of the 3-rd and 4-th 
moments change only a few %. Such a weak density dependence is easily understood 
in the chiral limit. If we take the zero current mass m = 0, the resulting pion mass 
becomes zero = = due to the Goldstone boson nature. Inserting = into 
(|10|) and p^), functions (7,^99 and I3 are rewritten as, 



] Jp = 0) = --r^iP hip ) / . 1 5- 



/3(/ = 0)-i/ ^ i^Jl-^ . (26) 

Namely, (7"^^^ and /a are reduced to the same function in the chiral limit. Using 
this result with eqs. (ET]) and (P3|), we find the 3-rd and 4-th moments become just 
numerical numbers, 

^val _ i + ^ _ 1 (27) 

^ 4 12 3 ^ ' 

j^val _ 1 + 1 _ 1 (28) 

and thus density dependence of both moments vanishes. Therefore, the valence quark 
distribution function in the chiral limit (zero mass pion) is independent of the matter 
density. In addition, from the results (|T^ |19|, |2^, [28| ), we can write an analytical 



expression of the quark distribution for the zero mass pion, xq{x) = x (0 < x < 1) 
at the model scale. 

In the case of SET2, the medium effect on the distribution function is drastic as 
clearly seen in Fig. 4 and 5. For x > 0.6, the result at the nuclear matter density 
shows a substantial reduction. In this value of the 3-rd moment decreases 

about 10%, and 20% for the 4-th moment. 

Finally, we discuss a connection between the nuclear pion enhancement of the sea 
quark distribution and the medium modification of the pion structure function. From 
eq. (0) calculated by the random phase approximation, it can be shown that the pion 
momentum distribution in the nucleus /^Ar(l/) is peaked at small y; y < 0.3 1^ |10 . 
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Thus, the integral of eq. ([l|) is dominated by the small y contribution. Let us recall 
that argument of q-n{x.„ = x/y) is proportional to 1/y. Hence, the integral (|l]) is 
sensitive to the large behavior of the quark distribution function||24||. We have 
shown that the quark distribution function at the large x^^ decreases as the density 
increases. Hence, we expect that the reduction of the pion structure function at 
the large Bjorken-x in the medium compensates the enhancement of the sea quark 
distribution due to the pion excess in the nucleus. In fact, it can be shown that, for 
X > 0.1, the enhancement of the sea in the nucleus almost disappears by using the 
density dependent pion structure function of SET2. On the other hand, the nuclear 
pion enhancement of the sea quark still remains in the case of SETl. More detailed 
discussions will be given in the forthcoming paper [E^ . 



In conclusion, we have studied the pion structure function in the nuclear medium 
within the NJL model as an effective theory of QCD. We have calculated the mo- 
ments of the pion structure function based on OPE to separate the soft part of the 
matrix elements and to maintain the gauge invariance. The result for the free pion 
shows a reasonable agreement with the experiment. We have shown that the quark 
distribution at the large x region decreases in the nuclear medium p = pq. Magni- 
tude of the modification depends on treatments of the medium effect. We have found 
that the medium modification of the pion structure function is small within the NJL 
model around the nuclear matter density. Instead, if we use the Brown-Rho scaling 
relation, the quark distribution function shows substantial reduction for the large x. 
We note that the quark distribution function in the high density limit approaches the 
non-relativistic one, q{x) = 6{x — 1/2) (at the model scale). It is also interesting to 
study the medium modification of other meson structure functions, which is under 



consideration 24 



We have also discussed the nuclear pion enhancement of the sea quarks in the 
nucleus in terms of the density dependent pion structure function. The reduction of 
the pion structure function at the large x may lead a suppression of the sea quark 
distribution due to the pion cloud. Such a tendency seems to be consistent with 
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the available experiments, where the sea quark enhancement in the nucleus is not 
observed. 
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Figure Captions 



Fig. 1 

The nucleon-hole and delta-hole excitations contribute to the deep inelastic lepton 
nucleus scattering. Solid line denotes the nucleon, and the double one the delta. The 
pion and virtual photon are depicted by the dashed and wavy curves, respectively. 

Fig. 2 

Graphical representation of the operator product expansion. The quark and the pion 
are depicted by the solid and dashed curves, respectively. Wavy curve denotes the 
virtual photon. Left hand diagram shows the forward scattering amplitude of the 
pion in the NJL model. 

Fig. 3 

The valence quark momentum distribution function at = 8GeV^ as a function 
of the Bjorken x. The solid curve denotes theoretical calculation for the free pion 
(p = 0), and the dashed curve the experimental fitp3|. 

Fig. 4 

The valence quark distribution function of the pion in the nuclear medium p = po ai 
g2 = sGeVl The result of SETl is depicted by the solid curve, and SET2 by the 
dashed curve. See text for detail. 

Fig. 5 

Ratio of the quark distribution functions at the normal nuclear matter density [p = 
Po) to one of the free pion (p = 0); 5'(a;)p=po/g(x)p=o- Results of SETl and SET2 are 
depicted by solid and dashed curves, respectively. 
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